Competing mechanisms of chiral symmetry breaking 
in a generalized Gross-Neveu model 
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Chiral symmetry of the 2-dimensional chiral Gross-Neveu model is broken explicitly by a bare 
mass term as well as a splitting of scalar and pseudo-scalar coupling constants. The vacuum and light 
hadrons — mesons and baryons which become massless in the chiral limit — are explored analytically 
in leading order of the derivative expansion by means of a double sine-Gordon equation. Depending 
on the parameters, this model features new phenomena as compared to previously investigated 4- 
fermion models: spontaneous breaking of parity, a non-trivial chiral vacuum angle, twisted kink-likc 
baryons whose baryon number reflects the vacuum angle, crystals with alternating baryons, and 
appearance of a false vacuum. 
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Consider the Lagrangian density of N species of mas- 
sive Dirac fcrmions in 1+1 dimensions with attractive, 
U(iV) invariant scalar and pseudoscalar interactions, 
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Flavor indices are suppressed (tpip = X)fc=i V'fcV'fc etc.) 
and the large N limit will be assumed. This 3-parameter 
field theoretic model generalizes the (massive) chiral 
Gross-Neveu (GN) model [l| to two different coupling 
constants. Its massless 2-parameter version is related 
to the early work of Klimenko Q and has only recently 
been investigated comprehensively Jj|. Our main moti- 
vation for considering the Lagrangian ([1]) is to study the 
competition of two different mechanisms of explicit chiral 
symmetry breaking, both of which are well understood in 
isolation. The first one is kincmatical and familiar from 
gauge theories — the bare fermion mass. The second one 
is dynamical — breaking chiral symmetry through the 
interaction term while preserving parity. This seems to 
have no analogue in pure gauge theories. In the present 
work we do not attempt a complete solution of the model 
([T]) which would require extensive numerical computa- 
tions. To get a first overview of its physics content, we 
focus on the vicinity of the chiral limit at zero tempera- 
ture, where everything can be done in closed analytical 
form. 

Following 't Hooft [3, the large TV limit is implemented 
by letting N — > oo while keeping Ng 2 and NG 2 constant. 
As is well known, this justifies the use of semiclassical 
methods [J [U • Thereby the Euler-Lagrange equation of 
the Lagrangian ([1]) gets converted into the Dirac-Hartree- 
Fock equation, 

(i 7 %-5-i 75 P)V = 0, (2) 
where the scalar and pseudo-scalar mean fields are related 



to condensates (ground state expectation values) through 



P = -G 2 (^i 75 V)- 
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Further simplifications arise if we concentrate on static 
problems in the vicinity of the chiral limit, where the 
potentials are slowly varying in space. This allows us 
to invoke a systematic expansion in derivatives of S and 
P without assuming that the potentials are weak 0, 0] . 
As a result, we arrive at an effective bosonic field theory 
in which the Hartree-Fock potentials appear as complex 
scalar field (written here in polar coordinates), 



S - iP = pe l9 



(4) 



Note that this method can only handle full occupation 
of single particle levels at present. It was pioneered in 
Ref. [8[ and applied systematically to two variants of the 
Lagrangian ([1}, the massive chiral GN model (g 2 = G 2 , 
Ref. Ml) and the massless generalized GN model (too = 0, 
Ref. Q ) . Since the form of the Hartree-Fock equation is 
the same in all of these cases, the problem at hand differs 
from previous ones only through the form of the double 
counting correction to the energy density, 
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An irrelevant term 



has been dropped. Regulariza- 
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tion and renormalization require only a straightforward 
extension of previous works. We replace the 3 bare pa- 
rameters (to ,(7 2 ,G 2 ) by physical parameters (^i,^2tV) 
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The In A dependence is mandatory to ensure that the 
ultraviolet divergence in the sum over single particle en- 
ergies is cancelled by the double counting correction. In 
the last line of Eq. (j6]) , we avoid the use of the standard 
confinement parameter Q 
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at this stage. This is done in order not to mix the param- 
eters of the model with dynamical quantities, which may 
lead to confusion in the present 3-parameter model. Re- 
stricting ourselves to the leading order of the derivative 
expansion, we assume furthermore that the radius p is 
fixed at the dynamical fcrmion mass and that the chiral 
angle field 9 is slowly varying. These assumptions can be 
justified by looking at higher order terms of the deriva- 
tive expansion, but they also have a very simple physical 
basis: Close to the chiral limit, the would-be Goldstonc 
field 9 (the "pion" field) is the only one which can be 
modulated at low cost of energy @. The renormalized 
ground state energy density (including the vacuum con- 
tribution) corresponding to Lagrangian (jTJ) then reads 
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Fcrmion number is given by the winding number of the 
chiral field 
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All we have to do is to minimize the energy J dx£ clas- 
sically. As a result, we will get information on the vac- 
uum and its symmetries, as well as on light mesons and 
baryons in the vicinity of the chiral limit. For a homoge- 
neous vacuum, the truncated derivative expansion is ex- 
act since the condensates are spatially constant. Hence 
our results for the vacuum may be taken as the large TV 
limit without any further approximation. Light hadrons 
are those which become masslcss in the chiral limit. Here 
the derivative expansion can be viewed as a kind of chiral 
perturbation theory, reliable close to the chiral limit. The 
expression for the pion mass for example is of th e ty pe 
of the Gell-Mann, Oakes, Renner (GOR) relation [1J| in 
the real world. The fact that baryons emerge from a non- 
linear theory for the pion field with the baryon number 
as topological winding number is of course reminiscent of 
the Skyrme model in 3+1 dimensions (ill. Il2j. 

We first determine the vacuum. To this end, we mini- 
mize 2tt£ , Eq. ([5]). with respect to (cc-independent) p and 
9 — the dynamical fermion mass and chiral vacuum angle. 
This yields the transcendental equations 
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FIG. 1: Qualitative shapes of effective potentials, Eq. 
in the (£,7) half plane. Each inserted plot shows £(9) in the 
interval [— so that the endpoints have to be identified. 
The origin (£ = 0, 7 = 0) is the chirally symmetric point 
where the effective potential vanishes. When crossing the 
critical lines 7 = ±£, the number of extrema changes. 



Their solution requires a case differentiation. To under- 
stand qualitatively what to expect, let us temporarily 
choose units such that the dynamical fermion mass is 1 
(p = 1) and focus on the ^-dependent part of the vacuum 
energy density, 

2nx£{6) = -27cos0- -£cos(20), £ = 6 - £1. (11) 

We have used the fact that the distinction between 77 and 
the confinement parameter 7 disappears in these units, 
cf. Eq. (|7||. Note also that depending on the sign of £, 
either the scalar coupling (for £ > 0) or the pseudoscalar 
coupling (for £ < 0) dominates. 

A survey of the ^-dependence of this effective potential 
in the (£, 7) half plane (7 > 0) reveals a rich landscape 
(see Fig. [1]): At the origin (£ = 0, 7 = 0), the potential 
is identically zero (not shown in Fig. 1) and the vacuum 
infinitely degenerate. This is the U(l) chirally symmetric 
point. Along the 7 axis there is a minimum at 9 = and 
a maximum at 9 = tt — the massive chiral GN model. 
Along the £ axis, there are two degenerate minima sepa- 
rated by two degenerate maxima — the massless general- 
ized GN model. As discussed in Ref. [|| , the minima can 
be identified with and tt for both £ > and £ < by 
means of a global chiral rotation, so that the positive and 
negative £ half-axes are in fact equivalent. What happens 
in the parameter region away from the 7- and £-axes de- 
pends on the sign of £. If £ > 0, the quadratic maximum 
becomes a quartic maximum at 7 = £; for larger values 
of £, a false vacuum develops at 9 = tt. In the limit 
7 — > the two minima become degenerate. If £ < on 
the other hand, the quadratic minimum becomes quartic 
when crossing the critical line 7 = — £. This is indicative 
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of a pitchfork bifurcation with two symmetric, degenerate 
minima present for 7 < — £. A non-trivial vacuum an- 
gle signals a non- vanishing pscudoscalar condensate and 
hence a breakdown of parity. This breakdown of parity 
is spontaneous, but induced by the explicit breaking of 
chiral symmetry. 

With this overall picture in mind, we return to 
Eqs. (fTUj) and solve them in two distinct cases: 

• Unbroken parity (phase I) 
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Broken parity (phase II) 
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In Eqs. Q12[) we have introduced the Lambert W function 
with the defining property 



x = W{x)e w{x) . 



(14) 



The vacuum energy in the parity broken phase II is lower 
than in the symmetric phase I. However, phase II only 
exists if # vac is real or, equivalently, 



6-6> W^ 1 )- 



(15) 



The next steps can be further simplified as follows. 
After minimization and determining the phase on the 
basis of Eq. (p~5]) . we normalize the radius of the chiral 
circle (the physical fermion mass) to 1 by a choice of 
units, 



p = p v 



1. 



(16) 



Then r\ may be identified with the confinement parameter 
Q familiar from the standard massive GN models, 



V = PI -> 7- 
In phase I, the condition p — 1 implies 

6 = 7- 

The vacuum energy density becomes 

4tt 2tt' 



(17) 
(18) 
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in agreement with the standard massive GN models. The 
^-dependent part of the energy density will be needed for 



the analysis of light mesons and baryons; in phase I it is 
given by 

1-k£1 = X -{B'f - 2 7 cos 9 - ^(6 - 7 )cos(2fl). (20) 
In phase II, the condition p = 1 implies 

6 = 0, (21) 
whereas the vacuum energy density assumes the form 
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In this phase, the ^-dependent part of the energy density 
reads 



27r£ e // = \{9'f - 2 7 cos 9 + i& cos(20). (23) 



Eqs. (|20I23|) can be treated simultaneously by setting 

2ir£ g = ^{9') 2 - 2 7 cos(9 - ^cos(26>) (24) 

with the definition 

re c _ / 6 - 7 (phase I, £ > -7) , . 

"> 2 ~ « - \ -fc (phase 17, C < -7) ( J 

We have traded the original bare parameters g 2 ,G 2 ,r] 
against two dimensionless parameters 7, £ and one scale, 
the dynamical fermion mass p = 1. The notation is cho- 
sen so as to agree with previous results for the massive 
chiral GN model Q for £ = and the massless general- 
ized GN model for 7 = 0. 

Next consider the light meson mass in both phases. 
Expanding expression (|24| around the vacuum angle 9 vac 
to 2nd order in •& = 9 — 9 V&C , we can simply read off the 
pion mass as follows: 

• Phase I (£ > -7) 







2tt£ ps -(i9') 2 + (7 + 6^ 2 + const. 
ml = 4(7 + 6 



(26) 



Phase II (£ < -7) 
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FIG. 2: Mechanical model illustrating vacuum structure, 
symmetries and meson masses of the generalized massive GN 
model (£ < 0) in the vicinity of the chiral limit. 



The last lines of Eqs. (|26I27[) may be regarded as the 
generalized GOR relations in our model. 

It is amusing that a well-known mechanical system is 
closely analogue to the present problem in the case £ < 0, 
cf. Fig. [3] A bead (mass m) is sliding without friction on 
a circular hoop (radius R) in a homogeneous gravitational 
field. The hoop rotates with constant angular velocity w 
around a vertical axis through its center. The Lagrangian 
reads 



L = \mB? ( 9 l + w' z sin^ 9 ) + mgR cos 9. 



(28) 



Denote the pendulum frequency by luq = \J g / R. At 
w = 0, there is a unique stable minimum at 9 = 0, ac- 
companied by small oscillations of frequency ojo- If one 
increases w, this minimum stays stable at first, but the 
frequency decreases like ^/ us^ — uS 1 until it vanishes at 
the critical value u = uq. At this point, two symmetric 
stable minima at 6 = iarccos (w 2 /w 2 ) develop, a text- 
book example of a pitchfork bifurcation [13|. Beyond 
this p oint, the frequency of small oscillations is replaced 
by \J u A — luq/lu. The gravitational field and the uni- 
form rotation are two distinct mechanisms of breaking 
the original SO(2) symmetry of the circle. The mapping 
of this mechanical problem onto our field theory model 
is obvious: U(l) chiral symmetry corresponds to the ro- 
tational symmetry of the circle, the bare mass plays the 
role of gravity, the difference in coupling constants cor- 
responds to the uniform rotation, the pion masses to the 
frequencies of small oscillations. We only have to identify 
7 = Wq/4, £ = — to 2 /A to map the two problems onto each 
other quantitatively. In principle, the regime £ > could 
also be modeled by assuming that the particle is charged 
and invoking an additional constant magnetic field, but 
in the absence of a phase transition this is less instructive. 

Let us now turn to baryons and baryon crystals. Here 
we need large amplitude solutions of the equation 



9" = 47sin6> + 2£sin2(9. 



(29) 



For small values of the parameters £, rj the kink-like soli- 
ton solutions of this equation are slowly varying so that 



the derivative expansion is applicable. The same is true 
for periodic soliton crystal solutions at sufficiently low 
density. However there is no restriction on the ratio £/7, 
so that the full phase structure shown in Fig. [T] is acces- 
sible in the vicinity of the chiral limit. Since Eq. (|2"9")l has 
no explicit x- dependence, it can be integrated once, 



1'f +47cos6» + £cos(26») = const. 



(30) 



The second integration is then carried out by separation 
of variables. 

The mechanical interpretation of the kinks is well- 
known: If we interpret x as time coordinate, Eq. (|29[) 
describes motion of a classical particle in a potential in- 
verted as compared to the potentials shown in Fig.[TJ The 
kink-like tunneling solutions between different vacua in 
field theory go over into classical paths joining two de- 
generate maxima in the mechanics case. In this classical 
mechanics interpretation, Eq. (|30p expresses conservation 
of the Hamilton function. As a matter of fact, Eq. (f2"9")) 
is nothing but the double sine-Gordon equation, a widely 
used generalization of the sine-Gordon equation to which 
it reduces if either 7 or £ vanishes. Its solutions can be 
found in the literature, see e.g. [l4j . so that we refrain 
from giving any details of the derivation. Since 9 is an 
angular variable, kinks do exist everywhere in the (£,7) 
half-plane. Inspection of the effective potentials of Fig. [1] 
then helps to understand the following results: 

For £ > —7 (phase I) there is only one kink solution 



'kink 



= — 2 arctan 



i y7sinh(2V£ + 72;) 



(31) 



We define the branch of the arctan such that 9 goes from 
to 2tt along the x axis. For £ < —7 (phase II) there are 
two different kinks depending on how one connects the 
minima along the chiral circle, 
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(32) 

Here, our choice of the branch of arctan is such that 9 
goes from —9 vajC to 9 vac for the small kink and from # vac 
to 2tt — 9 vac for the large kink. The baryon numbers 
B = Nf/N are 
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with 6* V ac from Eq. ([27]) with the + sign. The terms small 
and large refer to the chiral twist of the two kinks which 
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FIG. 3: Scalar and pseudoscalar potentials for the small kink 
baryon in the parity broken phase II with £ — — 1.37 as a func- 
tion of z = m-nX. The straight lines are the asymptotic values 
coinciding with the vacuum condensates. There are two de- 
generate vacua with equal scalar and opposite pseudoscalar 
condensates, related by a parity transformation. 
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FIG. 4: Same as Fig. [3] but for the large kink baryon. 



in turn is reflected in the baryon number. The baryon 
numbers of a small and a large kink add up to 1 simply 
because these kinks correspond to the 2 possibilities of 
travelling from one minimum to the other one along a 
circle. Eqs. (|3U33[I refer to kinks with positive baryon 
number. By changing the sign of 9, these can be con- 
verted into antikinks with opposite baryon number. 

In Figs. |3] and IU we illustrate the scalar and pseu- 
doscalar potentials for the small and large kinks in the 
parity broken phase II. To understand these graphs, we 
recall that the two vacua are characterized by the chiral 
angles ±# vac , Eq. (|27|) . The parity even, scalar vacuum 
condensate (cos(9 V ac) is the same in both vacua, the par- 
ity odd, pseudoscalar condensate (— sin(9 vac ) has oppo- 



FIG. 5: Same as Fig. [3] but for the kink baryon in the parity 
restored phase I and £ = — O.77. The pseudoscalar vacuum 
condensate vanishes. 



site sign. This is reflected in the asymptotic behavior of 
S and P for the kinks which connect these two vacua. To 
contrast this behavior with baryons in phase I (unbroken 
parity, £ > —7), we show in Fig. [5] the kink baryon from 
Eq. (|31|l where now both S and P are periodic. For the 
parameters chosen here, it resembles closely the standard 
sine-Gordon kink. 

Note the following limits: 

• Massive NJL model (7 > 0,£ = 0): There is a 
unique minimum at 9 = 0. We recover previous 
(sine-Gordon) results @, H| with the help of the 
identity 



9 = =f2 arctan 



1 



±4 arctan e 2 ^. (34) 



• Massless generalized GN model (7 = 0,£ > 0): 
There are 2 degenerate minima at 9 = 0, tt and 
correspondingly 2 kink baryons with baryon num- 
ber 1/2. The limit is singular (see Fig. [6}: As 
7 — > 0, the kink develops a plateau which becomes 
infinitely wide at 7 = 0. The kink decouples into 
2 half- kinks each carrying baryon number 1/2 0]. 
As one sees in Fig. [TJ this happens when the max- 
ima in the inverted potential become degenerate or, 
cquivalently, the false vacuum and the true vacuum 
in the original potential become equal. 

It is worth mentioning that there is yet another soli- 
tonic solution of some physics relevance: If one is inter- 
ested in the decay of the false vacuum, one has to consider 
tunneling through the barrier. This in turn is related to 
the kink-antikink which starts from the lower maximum, 
is reflected at the barrier and returns to the starting point 
(the bounce [Hj]). Since we are mainly interested in the 
vacuum and low-lying hadrons here, we do not go further 
into this problem. 
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FIG. 6: Behavior of #kmk as a function of z — m^x, for y/£ = 
KT 1 , 1(T 3 , 1(T 5 , 1(T 7 , with increasing width of the plateau 
for decreasing ratio 7/^. The plateau becomes infinitely wide 
in the limit 7 — > 0, leading to decoupled kinks with baryon 
number 1/2 in the massless model 0]. 



We now turn to a useful test of the consistency of our 
results, following Ref. [ll|. Consider the divergence of 
the axial current as obtained from the Euler-Lagrange 
equations for the Lagrangian ([!]), 



« = 2#y^ [m - (g 2 - G 2 )^] 



(35) 



The right-hand side exhibits the 2 sources of chiral sym- 
metry breaking, the bare fermion mass and the splitting 
of the coupling constants. The self-consistency condi- 
tions ([3]) and the renormalization scheme ([6]) can be used 
to rewrite Eq. (f3"5"j) as 



2NP 



[v - (6 - 6)5] 



(36) 



or, in units p = 1 and with the notation of Eqs. (|17I25[) 

2NP 



(37) 



Taking the expectation value of this equation in a time- 
independent state and remembering that j\ = j in 1+1 
dimensions, we arrive at the following expression for the 
fermion density, 



f(x) = / dx'P(x')[ 1 + ^S(x')} 

n J-00 

and, after another integration, the sum rule 
2N f°° 

N f = / dxxP(x) [7 + £S(ar)l . 
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The last equation in particular provides us with a non- 
trivial way of testing the baryon potentials. By inserting 



S = cos 9 and P = — sin 9 into the sum rule with 9 from 
Eqs. (|31l32p . we indeed reproduce the baryon numbers 
([33]l . Notice also that the expectation value of Eq. (|36|) 
for the divergence of the axial current, 



2NP(x) 



(40) 



reduces to the double sine-Gordon equation, Eq. (f2U)) , if 
we insert 



(41) 



and express S, P in terms of the chiral angle 9. This 
points to an alternative derivation of the basic equation 
(f29| which would not even require the derivative expan- 
sion, at least to leading order considered here. 

It is straightforward to compute the baryon masses by 
integrating the energy density and subtracting the vac- 
uum contribution, 



2ttM = 
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and the same results for the corresponding antikinks. 
These expressions are of course known from studies of 
the classical double sine-Gordon equation. 

Finally, consider baryonic matter at low density. The 
pertinent solutions of the double sine- Gordon equation 
are kink crystals which can be evaluated analytically in 
terms of Jacobi elliptic functions. Since we work only to 
lowest order of the derivative expansion in the present 
study, we bypass the complicated exact solution by sim- 
ply gluing together kink solutions. This is adequate in 
the low density limit. In the parity preserving phase I, 
the basic building block is #kink, Eq. (|3"Tj) . Let us denote 
the separation between two kinks (i.e., the lattice con- 
stant) by d, so that the baryon density is ps = 1/d. A 
dilute periodic array of kinks is then well approximated 
by 



crystal 



— Pkink 



0=- 



- nd) + 2-Kn for x € [nd — d/2,nd + d/2]. 

(44) 

For sufficiently large d this yields a smooth staircase 
curve which solves the double sine Gordon equation ex- 
actly except at the gluing points x = (n + l/2)d. There 
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FIG. 7: Unit cell of soliton crystal in the parity broken phase 
for 7 = 0.2, £ = —0.26. Small and a large kinks must alternate 
in the crystal due to their different asymptotics. The baryon 
numbers of the 2 constituents in the unit cell add up to 1. 



the error can be made arbitrarily small for large d. The 
energy density in the dilute limit is just MkmkAB with 
the kink mass from Eq. (|43|) . In phase II, we have to 
proceed slightly differently. Obviously one can only glue 
together the small and large kinks in an alternating way, 
see Figs. 13141 We therefore first construct a unit cell of 
the crystal by joining one small and one large kink, 



''kink 



Or) 



#small(# 
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-d/4) 
d/A) 



for 
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-d/2 <x<0 , 
0<x<d/2 



This carries baryon number 1 and is periodic modulo 2ir, 
so that the unit cells can now be assembled into a crystal 
in the same way as in phase I, Eq. (j33J), 

^c/ystai = ^kink(a^ — nd) + 2im for x £ [nd~ d/2, nd + d/2}. 

(46) 

The energy density in the low density limit of phase II 
becomes 



An example for a unit cell is shown in Fig. [7] with the 
same ratio £/7 and hence the same shape of the small 
and large kinks as in Figs. 13141 

Summarizing, we have investigated a 3-parameter gen- 
eralization of the U(l) chirally symmetric GN model. 
The two dimensionless parameters 7 and £ stem from two 
different mechanisms of breaking chiral symmetry explic- 
itly, the bare mass term and the difference between scalar 
and pscudoscalar couplings. Close to the chiral limit, the 
leading order derivative expansion has revealed the fol- 
lowing scenario. If the scalar coupling dominates, we find 
in general a unique vacuum with scalar condensate, light 
pions and kink-like baryons with baryon number 1. In the 
region £ > 7 a false vacuum shows up in the form of a 
second local minimum. If the pscudoscalar coupling dom- 
inates, at first nothing changes. Starting from a critical 
strength of the coupling (£ < —7), two symmetric min- 
ima appear together with scalar and pscudoscalar con- 
densates; parity is spontaneously broken. The mechani- 
cal model of a particle on a rotating circle in the gravita- 
tional field illustrates nicely the concomitant pitchfork bi- 
furcation. The two ways of connecting two minima along 
the chiral circle are reflected in two baryons whose baryon 
numbers add up to 1. These chirally twisted baryons are 
mathematically well known from studies of the double 
sine-Gordon equation and quite different from another 
typ e of twisted bound state specific for the chiral limit 
EH]- I n our case i the baryons are stabilized by topol- 
ogy. Shei's bound state is stabilized by partially filling 
the valence level and does not carry baryon number as a 
result of a cancellation with induced fermion number [l6| . 
In many respects the limits 7 — » and £ — > are atyp- 
ical so that previously explored 2-paramctcr versions of 
the present model cannot convey the full picture of chiral 
symmetry breaking in 4-fermion models. In view of the 
rich structure of the 3-parameter model, it seems worth- 
while to pursue its study, in particular to explore the fate 
of the symmetries at finite temperature and chemical po- 
tential. 



£ = (A/ sm all + Marge )pB 



(47) 



where the sum of the kink masses from Eqs. (|43[) can be 
simplified to 



small 
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